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CROSS-SPECTRAL  PROPERTIES  OF  SOME  COMMON  WAVEFORMS 
IN  THE  PRESENCE  OF  UNCORRELATED  NOISE 

INTRODUCTION 


An  extensive  amount  of  data  analyst  a had  been  performed  in  1974  to 
determine  towed  array  self-noiae  mechanisms.1*2*3  Reduction  of  the  data, 
with  particular  emphasis  placed  on  determining  the  magnitude  coherence  and 
phase,  yielded  the  cross-spectral  properties.  Since  magnitude  coherence  can 
be  a measure  of  aignal-to-noise  ratio*  and  in  conjunction  with  the  phase  is  a 
measure  of  wave  speed,  it  was  possible  to  determine  some  of  the  dominant  self- 
noise  mechanisms.  Interpretation  of  the  data,  however,  is  not  always  clear 
cut.  In  the  presence  of  uncorrelated  background  noise,  various  mechanisms  can 
combine  in  many  different  ways  to  produce  complicated  results.  To  achieve  a 
better  understanding  of  the  self-noise  data,  it  was  found  useful  to  investi- 
gate the  cross-spectral  properties  of  some  common  waveforms  in  the  presence  of 
uncorrelated  noise.  It  is  the  purpose  of  this  report  to  present  the  results  of 
this  investigation,  with  the  hope  that  it  will  prove  useful  to  others  analyzing 
these  types  of  data. 

THEORETICAL  BACKGROUND 


This  section  provides  the  theoretical  background  for  determining 
magnitude  coherence  and  phase.  For  a single  traveling  wave  in  the  presence  of 
uncorrelated  noise,  Carter*  showed  that 


y (f)l  * (1) 

T12U;  X(f)  + N(f ) * u# 

based  on  the  model  in  figure  1,  where  | F_ 2 ls  the  magnitude  coherence 

between  observation  points  1 and  2.  Uncorrelated  noise  sources  ni(t)  and  n2(t) 
have  the  respective  autospectra  N^(f)  and  No(f),  where  N^(f)  ■ N2(f)  ■ N(f). 
X(f)  is  the  autospectrum  of  x(t),  which  is  100  percent  coherent,  since  it  is 
the  same  at  both  observation  points.  (A  computer  run  experiment  using  actual 
at-sea  data  verified  equation  (1)  and  is  further  discussed  in  a NUSC  memoran- 
dum.^) 

n|(0 


x(0 


DELAY  n2(f) 


SENSOR  1 


SENSOR  2 


Figure  1.  Model  for  Single  Traveling  Wave 
Plus  Uncorrelated  Noise 


For  two  independent  traveling  waves  in  the  presence  of  uncorrelated  noise, 
Gardner 5 has  shown  that 
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Yl2<0  " 


rx<f> 

f*n<»  x22{f>] 

e10x(f)  + ry(f) 

YU<f)  Y22(f) 

>*  e 

Jxn(f)  + Yn<Oj 

X22(f)  + Y22(f j 

V 

(2) 

where  Xj^f),  X22(f),  V1,(f),  and  Y22(f)  are  the  autospectra  for  the  two  inde- 
pendent signals  fx(t)  ana  y(t))  at  observation  points  1 and  2.  Here  uncorre- 
lated noise  is  included  as  part  of  the  signal*  The  cross-spectral  densities 
are  rx(f)[X1}(f)  X^Cf)]1*  e1**^)  and  ry(f)[Yn(f)  Y^f)]**  e1*^) 

(i.e.,  X12(f5  and  Y12(f),  respectively)  and  the  phases  are  0x(f)  ant 
Equation  (2)  can  be  written  as 


and  0y(f). 


Y12(f) 


X(f)e10x(f)  + Y(f)ei0y(f) 
X(f ) + Y(f)  4 N(f  ) 


(3) 


as  shown  by  Markowitz0,  for  the  model  In  figure  2 . 

nx1(n 


*(t) 


y(t)  ■ 


TIME 

DELAY 


TIME 

DELAY 


ny)(0 


ny2(0 


T 

nx2(f> 


SENSOR 

1 


SENSOR 

2 


Figure  2.  Model  for  Two  Traveling  Waves  Plus  Uncorrelated  Noise 


In  the  figure  n^t),  n„2(t),  nyi(t),  and  ny2(t)  are  uncorrelated  noise 
sources  with  autospectra  Nxj^(f)  - Nx2(f)  “ Nx(f)  and  Nyi(f)  “ Ny2(f)  - Ny(f). 
The  signals  x(t)  and  y(t)  are  each  100  percent  coherent  in  the  model  and  have 
the  respective  autospectra  X(f)  and  Y(f).  Also,  Xj^(f)  ■ X22(f)  "X(f)  + Nx(f), 
Yn(f)  - Y22<f>  " Y(£)  + Ny(f),  and  Nx(f)  + Ny(f)  - N(f).  We  note  that  (flx  « 
u>d/cx)  and  0y  “ (ud/cy),  where  d is  the  spatial  separation  between  sensors  1 
and  2 and  cx,  cy  are  the  respective  x and  y velocities. 

A derivation  of  equation  (3)  is  given  in  the  appendix.  The  important 
assumptions  made  there  are  that  (1)  waves  are  nondispersive;  (2)  no  attenua- 
tion between  observation  points  1 and  2;  (3)  x(t)  and  y(t)  signals  in  equation 
(3)  are  independent  and,  therefore,  uncorrelated  with  respect  to  each  other; 
and  (4)  the  estimate  of  the  initial  phases  associated  with  each  wave  are 
uniformly  distributed  from  -it  to  +ir. 
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The  magnitude  coherence  and  phase  obtained  from  equation  (3)  are 

|x(f)2  + Y(f)2  + 2X(f  )Y (f ) coa  ^(f)  - 0y(f)J  | **  {4) 


y12(£)  « 


and 


X(f ) + Y(f ) + N(f) 


jX(f)  sin  0x(f)  + Y(f ) sin  0y(f) 

tan  (X(f)  cos  0 (f)  + Y(f)  cos  0y(f) 

x y 


(5) 


The  significance  of  the  third  assumption,  regarding  the  statistical  inde- 
pendence of  x(t)  and  y(t),  can  readily  be  seen  in  the  special  case  of  standing 
waves.  Consider  x(t)  and  y(t)  to  be  of  unity  magnitude  and  traveling  in  oppo- 


site directions  at  frequency  ui0  ■ 2irf 0 to  form  the  standing  wave  Sj  ” cos 
k xje*wot  at  observation  point  1 and  S2  m cos  k0x2eiwot  at  observation  po 


One  can  then  form  the  cross  correlation  to  obtain 

T/2 


point  2. 


R12(T)  “ \f COS  koXlCOS  koX2 


-i2nf  t i2irf  (t+t) 
s o e o dt 


-T/-2 

lim  T -*■  « 

cos  k x.  cos  k x,ei2,lfoT. 
o 1 o l 


The  cross-spectral  density  becomes 

OO 

n tc\  C 1 1 i2irf  t i2irfx, 

G10  (f)  ■ /cos  k x-  cos  k x0e  o e dx 
J 0 1 O L 


cos  k x,  cos  k x„  2tt6  (2itf-2irf 


o 1 


o 2 


cos  k x,  cos  k x06 
o 1 o 2 


K) 


Normalizing  via  C^(f)  and  G^ (f)  yields 


Y12(f) 


cos(u)  x,/c)  cos  (u  x./c) 
o 1 0 L 

|cos(w  Xj/c)|  jcos(u>ox2/c)| 


(6) 


0 in  phase 

0 - 

it  out  of  phase 


(7) 


where  c is  the  wave  speed  (to  /k  ) . (Uncorrelated  noise  sources  are  ignored  tem- 

“OO 

porarily,  since,  for  the  purpose  of  demonstrating  the  significance  of  the  third 
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assumption,  they  represent  an  unnecessary  complication.)  These  results  (equa- 
tions (6)  and  (7))  differ  substantially  from  the  previous  results  (equations  (4) 
and  (5)).  For  a standing  wave,  equations  (4)  and  (5)  show,  respectively, 

|y12(f)  | to  be  periodic  in  f and  0 to  be  spatially  homogeneous,  whereas  equa- 
tions (6)  and  (7)  show,  respectively,  |y^2^)I  “ 1 and  0 to  be  nonhomogeneous , 
depending  upon  actual  values  of  x.  and  X2>  T'.ie  difference  lies  in  the  fact  that 
in  equations  (6)  and  (7)  the  x ana  y traveling  wave  components  of  the  standing 
wave  are  not  statistically  independent,  but  in  equations  (4)  and  (5)  they  are. 
This  difference  is  shown  mathematically  in  the  appendix,  where  the  estimate  of 
the  crosd-spectral  properties  of  the  combined  signals  are  obtained  via  ensemble 
averaging.  In  each  computation  of  the  cross-spectral  density,  an  initial  phase 
is  associated  with  each  traveling  wave.  When  the  two  waves  are  statistically 
independent,  this  initial  phase  is  taken  to  be  random,  with  a uniform  distribu- 
tion from  -it  to  it.  When  the  ensemble  average  is  formed,  various  cross-product 
terms  drop  out  because  of  the  random  phase  relationship  and,  thereby,  equations 
(4)  and  (5)  are  obtained.  In  equations  (6)  and  (7),  the  initial  phase  is  assumed 
to  be  zero,  resulting  in  a deterministic  phase  relationship,  namely  zero;  conse- 
quently the  results  are  different. 


RESULTS 


GENERAL 


This  section  is  essentially  a graphic  representation  of  equations  (4)  and 
(5)  for  variations  of  the  pertinent  parameters.  These  parameters  are  relative 
wave  speed  determined  by  0X  and  0y,  relative  magnitude  determined  by  X(f)  and 
Y(f),  and  relative  uncorrelated  noise  determined  by  N(f).  For  simplicity,  the 
frequency  dependence  of  X(f),  Y(f),  and  N(f)  is  neglected;  i.e.,  the  autospectra 
of  x(t),  y(t),  and  n(t)  are  constant  with  respect  to  frequency,  and  equations 
(4)  and  (5)  become,  respectively. 


ri2(f> 


r 2 2 

'x  + Y + 2XY 
L 


cos(0  - 
x 


X + Y + N 


(8) 


and 


(X  sin  0 + Y sin  0 ) 

0 “ tan  -jti a,X  , „ 

(X  cos  0 + Y cos  0 ) 
x y 

where  X(f),  Y(f),  and  N(f)  are  now  simply  written  X,  Y, 


and  N. 


(9) 


Results  are  presented  as  a function  of  0X  for  various  values  of  0y/0x. 

In  actual  practice,  | y 12 1 and  0 are  obtained  as  functions  of  frequency.  The 
conversion  is  readily  made  since 

0 “ 2nfd/c  and  0 ■ 2ufd/c  , (10,11) 

x x y y 


where  d is  the  spatial  separation  between  sensors  1 and  2,  and  cx  and  cy  are  the 
respective  wave  velocities  for  x(t)  and  y(t). 

In  general,  the  coherence  function  will  reflect  a periodicity  with  respect 
to  frequency  because  of  the  term  [2XY  cos  2nf^x-  In  fact,  a single  period 
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with  respect  to  frequency,  denoted  A£,  is  exactly  cxc„/d(cv  - c„): 

* y y x 

Af  ■ cxcy/d(cy  " Cx)>  (12) 

This  periodicity  is  very  important  in  determining  self-noise  mechanisms,  since 
it  is  a reflection  of  wave  speed,  which  is  associated  with  the  self-noise  mech- 
anisms (l.e.,  5000  ft/sec  is  associated  with  acoustic  waves).  Of  course,  the 
fact  that  we  are  dealing  with  two  speeds  can  present  complicated  results.  When 
only  one  wave  speed  exists  (l.e.,  cx  ■ Cy  or  Y * 0)  Af  -*■  <»  and  | Y12I  *8  a constant. 

We  note  that  the  phase  is  independent  of  N.  However,  from  a computational 
point  of  view,  uncorrelated  noise  is  important  in  actual  practice,  since  the 
variance  for  the  computed  estimate  of  the  phase  depends  on  the  coherence  and, 
therefore,  on  N. 

SPECIFIC  CASES 

Discussion  here  is  limited  to  waveforms  resulting  from  the  combination  of  no 
more  than  two  traveling  waves.  Four  waveforms  are  investigated: 

1.  Single  Traveling  Wave  Plus  Uncorrelated  Noise 

Equations  (8)  and  (9),  respectively,  reduce  to 


|Y12<f)l  * X+N  and  0 " K 


(13,14) 


(see  figures  3 and  4).  We  note  that  Iyj^I  is  constant  and  this  indicates  that 
only  one  wave  speed  exists,  hence  only  one  self-noise  mechanism.  This  wave  speed 
is  c ■ 2nd/ (0/f),  where  (0/f)  is  the  slope  of  the  phase  plot. 

2.  Waves  With  0y  - -0X  Plus  Uncorrelated  Noise 


Equations  (8)  and  (9)  reduce  to 


lr12(f)  | 


f 2 2 

X + Y + 2XY  cos 


»„)' 


and 


tan 


X + Y + N 


[(X  - Y)  tan  0x] 

OTTy)  j • 


(15) 


(16) 


The  magnitude  coherence  is  periodic  in  0X,  with  a period  of  n radians  corre- 
sponding to  Af  ■ cx/2d.  The  speed  of  the  traveling  components  of  waves  with 
0 y ■ -0X  is,  therefore,  c - 2dAf. 

2a.  Equal  Magnitude  Traveling  Wave  Components  (Classical  Standing  Wave) 


Equations  (15)  and  (16)  become,  respectively, 


Iy12(oI 


— (2  + 2 cos  20x)Si 


2X 

2X  + N 


cos  0 

x 


2X  + N 


(17) 
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1.0 

■ 

y|T12j’0.5 

for  s,gnal  =t 

NOISE 
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Figure  3.  | Y12 1 

for  Single  Traveling  Wave 

Plus  Uncorrelated  Noise 

Figure  4.  0 for  Single  Traveling  Wave 
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and 


0 - 


3* 


0 £ Dx  < /2,  if  < tx  < /2  tan  0x  ia  positive) 


n/2  K #x  1 »•  3w/2  < 0X  < 2 it  (l.e.,  tan  0x  Is  negative) 


(18) 


The  essential  characteristics  of  the  classical  standing  wave, 
shown  in  figures  5-7,  are  as  follows: 


• Figure  5 shows  the  step  function  characteristic  of  the  phase 
where  0 contains  jumps  of  tt  radians.  This  is  superimposed  over 
the  traveling  wave  component  for  reference. 

• Figure  6 shows  the  periodic  characteristic  of  the  standing  wave 
for  2X/(2X  + S)  ■ 0.6.  The  nominal  values  chosen  were  X ■ Y ■ 3 
and  N ■ 4.  We  note  that  for  equal  magnitude  components  the 
minima  occur  at  | Y12 1 " 0 &nd  the  maxima  at  2X/(2X  + N).  The 
frequency  period  is  Af  - cx/2d  or,  with  respect  to  0x,  it  is  tt. 

• Figure  7 shows  that  as  the  noise  (N)  is  reduced,  the  peak  in- 
creases to  the  new  value  of  2X/(2X  + N),  which,  in  this  case, 
is  0.8.  The  nominal  values  chosen  were  X ■ Y - 2 and  N ■ 1. 

2b.  Unequal  Magnitude  Traveling  Wave  Components 

For  Y - 2X,  as  an  example,  equations  (15)  and  (16)  become,  re- 
spectively, 

X(5  + 4 cos  20  )**  X , , 

I Y12(f  > I “ 3X  + N ’ 3X  + N (1  + 8 cos  *x> 


and 


0 • tan-1  (1/3  tan  0x). 

These  : ..alts  are  shown  in  figures  8-10.  As  for  all  waves  with  0y  - -0X,  the 
periodicity  of  Iyj^I  with  respect  to  0X  is  n. 

The  minima  are  no  longer  zero,  as  previously  observed  for  X ■ Y, 
but  depend  on  the  noise  (N).  The  effect  of  N is  to  control  the  minimum  points, 
peak  points,  and  peak-to-valley  excursions.  These  values  decrease  as  N Increases. 

Characteristics  of  the  phase  plot  are  seen  in  figure  11  to  be 
between  the  standing  wave  step  function  (X  m Y)  case  and  the  straight  line 
single  traveling  wave  case,  which  again  is  shown  for  reference  in  the  figure. 
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Figure  6.  |y12I  for  Standing  Wave  (for  2X/ (2X  + N)  - 0.6) 


Figure  7.  | y12  | for  Standing  Wave  (for  2X/(2X  + N)  - 0.8) 
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Figure  11.  0 for  Two  Waves  With  0 ■ -0 

y x 


We  can  now  generalize.  For  waves  with  0 ■ -0  , the  range  of 

y ^ 

X - Y | I . X + Y 

X + Y + N -P12I-X  + Y + N * V ; 

and  if  we  assume  by  convention  X > Y,  then  the  mean  value  of  |y12I  wil1  be 

i Yi2 1 mean  “ x'Vy'+'n  X > Y-  (20) 

with  an  excursion  of  +Y/(X  + Y + N),  as  shown  in  figure  12. 

Therefore,  as  Y -►X,  for  a given  N,  the  minima  approach  zero,  the 
peak-to-valley  excursion  is  maximized,  and  | Y12 1 mean  reduced.  For  X >>  Y, 
the  results  approach  the  single  traveling  wave  case.  This  is  seen  in  figure  13 
for  X - 2Y,  X - 5Y,  X - ilY,  and  N - 2/3(X  + Y).  Thus,  as  (X/Y)  increases, 

| Y^2 1 becomes  more  nearly  constant. 

The  generalization  pertaining  to  0 is  simply  that  as  X/Y  1 
the  characteristic  step  function,  seen  in  figure  5,  is  observed,  and  as  X >>  Y 
the  results  approach  the  single  traveling  wave  case  shown  in  figure  4.  Figure 
11  presents  an  "in  between"  case,  where  X - 2Y.  Finally,  figure  14  compares 
X - 5Y  with  X - 2Y,  where  the  former  case  is  already  approaching  the  straight 
line  condition. 

3.  Two  Waves  Traveling  in  Same  Direction 

From  equation  (8)  we  note  that  the  magnitude  coherence  will  have  a 
periodicity  in  f dependent  upon  (0X  - 0y)  and  is  indistinguishable  from  other 
waveforms  with  similar  characteristics.  For  example,  a wave  -i^h  0y  **  - 0X 
having  traveling  wave  components  X8W  and  Yew  and  a phase  vel<  city  of  0xaw  is 
compared  with  the  sum  of  two  waves,  Xtw  and  Y^w»  traveling  in  the  same  direction 
with  phases  0xtw  an<*  flytw  Tbelr  magnitude  coherences  are  indistinguishable  if 


Iy 


12' 


is 
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0 3 • 6 * 


Figure  12.  Iy^I  ^or  Two  ^aves  With  0y  ■ -0X;  Maxima,  Minima,  and  Mean 


Figure  13.  |y1?|  for  Two  Waves  With  0y  - -0x;  Effect  of 
Relative  Wave  Magnitudes 


•x 


Figure  14.  0 for  Two  Waves  With  0y  “ -0X;  Effect  of 
Relative  Wave  Magnitudes 
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X8W  ■ X^,  Yw  ■ Yt  , and  20xgw  - (0xtw  “ ®ytw^*  This  situation  is  discussed 
further  on  p.  13. 

3a.  Equal  Magnitude  Components 

For  equal  magnitude  traveling  waves,  equation  (9)  becomes 

. fain  0 + sin  0 1 0+0 

• ■ tan  Lcos  0X  + cos  0yJ  - * (21) 

x y 

Thereby,  0 becomes  a linear  function  of  0X  for  given  values  of  (0y/0x)  and  is 
shown  In  figure  15. 


Figure  15.  0 for  Equal  Magnitude  Waves  Traveling  in 
Same  Direction;  Effect  of  0y/0x 


3b.  Unequal  Magnitude  Components 


In  general,  the  two  traveling  waves  will  have  unequal  magnitudes. 
In  such  a case,  certain  characteristics  of  0 can  be  observed  by  rewriting  equa- 
tion (9)  as 


0 - tan 


-1 


tan  0 


(22) 


In  the  above  equation,  we  assume  X to  be  the  dominant  component. 

(This  is  done  arbitrarily  for  convenience,  and  Y and  0y  can  easily  be  substituted 
for  X and  0X,  yielding  the  same  results.)  For  (X/Y)  » 1,  this  term  dominates 
and  0 0X,  regardless  of  the  va.'ue  of  (0y/0x).  For  (0y/0x)  < 1 and  X/Y  - 0(1), 

the  trigonometric  terms  introduce  periodic  variations  about  the  line  0 “ 0X.  At 
X/Y  ■ 1,  of  course,  equation  (21)  applies. 

In  actual  computations,  in  order  to  get  the  results  of  equation  (21), 
it  is  critical  that  X/Y  be  within  a fraction  of  a percent  from  1.  This  is  shown 
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in  figures  16  and  17.  Clearly,  for  (0  /0X)  ■ 1,  0 - 0X.  For  (0y/0x)  > 1 and 
X/Y  - 0(1),  again  periodic  variations  about  the  line  0 - 0X  are  introduced.  How- 
ever, computationally,  if  X/Y  gets  "too  close"  to  1,  jumps  of  it  or  2tt  may  occur. 
This  is  shown  in  figures  18  and  19.  Also,  the  periodicity  of  these  variations  is 
[l/(l-0y/0x]  2ir.  Thus,  in  figures  16  and  17,  the  periodicity  of  the  variation 
about  tne  line  0 - 0X  is  (10/9)2ir  and  4n,  respectively.  In  figures  18  and  19, 
the  periodicity  of  the  variation  about  the  line  0 ■ 0X  is  -2n  and  -u,  respec- 
tively. The  positive  or  negative  sign  determines  whether  the  variation  is  initi- 
ated above  or  below  the  line  0 “ 0X. 

4.  Two  Waves  Traveling  in  Opposite  Directions 

The  results  for  two  waves  traveling  in  opposite  directions  are  almost 
the  same  as  in  section  3.  (p.  10),  except  for  the  sign  change  in  0y  and  sin 
0y.  The  main  difference  is  in  the  periodicity  of  the  coherence  function  and  in 
the  periodicity  of  the  variations  about  the  line  0 - 0X.  As  mentioned  in  sub- 
section 3b.,  the  periodicity  of  these  variations  is  [ 1 / ( 1 - 0y/0x)]  2tt.  There- 
fore, for  (0y/0x)  “ -0.1  (see  figure  20)  and  (0y/0x)  - -2  (see  figure  21),  these 
periodicities  are  ( 10/11) 2n  and  2tt/3,  respectively.  For  comparison  purposes  one 
should  refer  to  figures  16  and  18.  There  the  periodicities  are  (10/9) 2tt  and  -2ir, 
respectively,  where  (0y/0x)  « + 0.1  and  + 2. 

ADDITIONAL  COMMENTS 

Unfortunately,  similar  and,  in  some  cases,  identical  phase  plots  arise  from 
completely  different  combinations  of  self-noise  mechanisms.  The  same  is  true  for 
coherence  plots.  These  ambiguities  are  usually  resolved  by  looking  at  both  the 
phase  and  coherence.  However,  care  must  be  exercised. 

As  a simple  example,  one  notes  that  for  equal  magnitude  traveling  waves 
equation  (21)  shows  0 = (0X  + 0y)/2.  Denoting  cases  1 and  2 by  their  respec- 
tive subscripts,  we  see  that  for  0X^  + 0y^  « 0X2  + 0y2  the  phase  plots  are  in- 
distinguishable. In  terms  of  wave  speed,  this  becomes 

di  p-  + -L\  - d2  (_L_  + _l\ 

\cxl  cy>)  \cx2  cy2/ 

Since  the  spacing  between  sensors  is  known  and  is  the  same  in  both  cases,  any 
combination  of  cxj,  cy^  and  cxo,  Cy2  that  satisfies  the  above  yields  the  same 
and,  therefore,  ambiguous  result.  Consider  an  investigator  who  observes  a 
straight  line  phase  plot  with  a slope  of  0.6°/Hz  for  a sensor  spacing  of  6 in. 

He  might  interpret  the  result  to  be  a single  wave  traveling  at  300  ft/sec,  where 
c - 360  d/(0/f).  Call  this  case  1.  However,  the  investigator  might  not  have 
noticed  electronic  noise  that  was  combining  with  a wave  equal  in  magnitude  to 
the  electronic  noise  and  traveling  at  150  ft/sec.  Call  this  case  2.  Then, 
cxl  “ cyi  “ 300  ft/sec,  while  cx2  ° 150  ft/sec  and  cy2  -+  °°.  Although  both  cases 

1 and  2 yield  the  same  phase  plot,  they  can  be  easily  distinguished  by  studying 
the  coherence  plots.  Case  1 shows  a constant  magnitude  coherence,  whereas  case 

2 reflects  a periodicity  with  respect  to  frequency. 

In  another  example,  a wave  with  0y  ■ -0X  is  compared  with  the  sum  of  two 
waves  traveling  in  the  same  direction.  The  wave  with  0y  = -0X  has  traveling 
wave  components  Xsw  and  Ysw  and  phase  0XSW.  The  two  waves  traveling  in  the 
same  direction  have  components  Xtw  and  Ytw,  with  phases  0xtw  ar>d  0ytw  It 
was  mentioned  in  section  3.  (p.  10)  that  the  magnitude  coherences  are  identical 
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(0. 


- 0 


—xsw  '•'xcw  "ycw^*  **  cxaw  “ 300  ft/sec, 
Cycw  “ 5000  ft/sec,  Che  results  are  not  identical  but 


if  - Xtw,  Ygy  - Ytw  and  20 
cxtw  “ 350  ft /sec,  and 

so  close  as  to  render  themindlstinguishable.  The  two  cases,  however,  are 
easily  distinguished  by  noting  the  phase.  A straight  line  drawn  through  the 
periodic  variations  will  determine  the  speed  of  the  dominant  wave  (in  this 
case  either  300  or  150  ft/sec). 


The  final  example  compares  two  waves,  denoted  Xg  and  Y8,  traveling  in  the 
same  direction  with  two  waves,  denoted  Xq  and  YQ,  traveling  in  opposite  direc- 
tions. Following  the  convention  adopted  in  subsection  3b.  (p.  12),  we  assume  X 
to  be  the  dominant  component.  Then,  if  X0  ■ Xg,  Y0  ■ Y8  and  if  cxo  “ cx8, 
the  only  observable  difference  in  both  the  phase  and  coherence  plots  would  be 
in  the  periodicity  with  respect  to  frequency  as  a result  of  the  fact  that  cyo  0 
Cy8.  The  observed  periodicity,  as  mentioned  previously,  is  [3/(l-0y/0x)}  2n 
with  respect  to  0jj.  Since,  in  actual  practice,  data  are  obtained  as  a function 
of  frequency  (not  0X,  which  J s what  the  investigator  is  seeking  to  determine), 
this  periodicity  is  Af  -fA/d\/^  - 3j)J  or,  denoting  (0y/0x)  - (cx/cy)  - a. 


Af  - c /d(l-a). 
x 


(23) 


In  order  for  the  periodicities  to  be  the  sane,  Afs  would  have  to  equal  Af0, 
which  would  require  that  a8  - a0,  since  d is  the  same  in  both  cases,  as  is  cx 
(by  definition).  This  is  impossible  because  <*8  is  always  positive  and  aQ  always 
negative.  It  is,  however,  possible  to  obtain  Af0  » -Af..  This  would  produce 
an  Identical  periodicity  in  the  coherence  plots  for  both  cases.  Also,  the  same 
periodicity  would  be  produced  in  the  phase  plots  for  both  cases,  but  with  a half- 
period  phase  shift.  This  is  all  shown  in  figure  22. 

It  is  hoped  that  the  results  obtained,  along  with  the  examples  presented, 
are  of  help  to  those  analyzing  these  types  of  data.  For  convenience  the  results 
are  summarized  in  table  1. 

In  actual  practice,  the  reduced  data  are  usually  more  complicated  owing  to 
(1)  the  presence  of  more  than  two  self-noise  mechanisms,  (2)  dispersive  waves, 

(3)  damping,  (4)  the  statistical  dependence  of  waves  as  a result  of  reflections, 
and  (5)  the  dependence  of  all  pertinent  parameters  on  frequency.  If  desirable, 
the  investigator  can  often  minimize  these  complications  by  careful  experiment 
design. 

The  investigator  should  also  look  at  the  time  domain  data,  in  particular 
the  cross  correlation,  to  shed  light  on  or  to  confirm  his  analysis.  Periodi- 
cities occurring  in  the  frequency  domain  might  be  seen  more  readily  via  a 
Cepstrum  analysis^  or  a Smoothed  Coherence  Transform  (SCOT) . *" 
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Figure  22.  Resolving  an  Ambiguity 
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Appendix 

DERIVATION  OF  MAGNITUDE  COHERENCE  AND  PHASE  FOR 
TWO  TRAVELING  WAVES  IN  THE  PRESENCE  OF  UNCORRELATED  NOISE 


TWo  traveling  waves  plus  random  noise  are  measured  at  observation  points 
1 and  2.*  At  u ■ u0  the  observed  signals  are,  respectively. 


X X 

lu  t ie  . la  (t-  — ) ie  lib  (t-  -i)  le 

„ o nl  o cv  x . „ o c y 

Sj^  - ^e  e +Xe  x e + Ye  Ye7 


and 


X X 

lb)  t 19n2  lb)  (t-  ) 10  lb)  (t-  ie 

S2  - n2e  0 e + Xe  ° C*  e X + Ye  0 "x  e y. 


(A— 1) 


(A-2) 


where  0„i,  0n2>  6x>and  0y  are  the  Initial  phases  associated  with  a particular 
observation,  and  nj  ■ n^(b)),  n£  ■ ^(w),  X ■ X(b>),  and  Y ■ Y(u)  are  understood. 
Also  note  X Is  a traveling  wave  as  distinguished  from  x the  spatial  coordinate. 

The  autocorrelation  at  point  1 is 

G:  2 ib)  t T(-2  ib>  t -lb)  A)  Ke-e  ,) 

Bj  e dt  + J n^Xe  e x e dt 

2 -T/2 


11m  T -*■  ® 
T/2 


T/2 

ib)  T -Ib)_x,/Cv  1(0  -9  , ) f - lb)  X 

y e y nl  dt  + / X2e  ° dt 


/iU)  T -J.O)  ~ . / 

V o O 1 
n^Ye  e 


-T/2 


/ 

-T/2 


T/2 

/x»i 

-T/2 


T/2 


T/2 


ib)  T lb)  x /c  1 <0  _ —0  ) /•  lb)  T lb) 

o o 1 x nix..  . / w„  o o 

e e e dt  + / XYe  e 


/lb)  T lb)  X, /c  1(0  ,-0  ) *'/•'  lb)  1 

YV  ° e ° 1 ■>  e 11  > d,t  / YXe  ° 

-T/2 
T/2  . 

r 7 iu)«T 
+ / Y2  e dt 


f 

-T/2 

T/2 

/ 

-T/2 


(H) 


i(0  -0  ) dt 
y x 


^ 

o\c„  c I 
! \y  x/e 


1(0  -0  ) 
x y 


dt 


£ 


(A-3) 


*The  approach  taken  in  this  appendix  was  suggested  by  Dr.  A.  H.  Nuttall,  of  NUSC. 
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The  corresponding  autospectrum  Is 


G. , (u)  - (N  2 + X2  + Y2)  2*6(w-w  ) 

11  1 o 

[-iw  x,/c  i(0  -6  ) iw  x,/c  i(0  .-0  )1 

e o 1 x e X nl  + e o 1 X e nl  x'J 

>iY[' 


2tt6(w-w  ) 
o 


-i<0  X, /c  i(0  — 0 _ ) iw  x,/c  i(0  ,-0 
+ n,Y  | e 0 1 y e Y «1  + £ o 1 y £ ' nl  y 


I 

fiw  1^1  - h\  i(0  -0  ) iw  i(0  -e  )1 

+ XY  [e  °\\  Cy)e  y x +e  o^Cy  ge  X yj 


ir6(w-w  ) 
o 


2it6(w-w  ).  (A-4) 

o 


Taking  0 , 0 0 yand  0 as  uniformly  distributed  random  variables  and 

ensemble  averaging^  we  ge? 


G. , (w)  - (n  2 + X2  + Y2)2n6(w-w  ). 
11  1 o 


(A-5) 


All  cross  terms  involving  differences  between  Initial  phase  angles  go  to 
zero.  This  can  be  seen  as  follows: 


i(e. -0  ) ?.  1.  Keu-ej 


b a' 


f f e b “ P(0b)P(0a)dV0a,  <A'6> 


-IT  -IT 


where  P(0,  ,0  ) - P(0,  )P(0  ) for  statistical  independence. 
d a d a 


Taking  P(0a)  - ^ 


-tt  < 0 < tt  and  P(0.  ) ■ — 

a b 2ir 


IT  < 0 < -TT 

a 


-TT  < 0.  < TT 

— D — 


7T  < 0 < — TT 

b 


we  get 


i(e. -o  ) 

d a 


7T  TT 


hff 


i(0.-e  ) 

e d0.d0  - 0. 

b a 


(A-7) 


Therefore,  P(0£)  or  P(8b)  must  be  uniformly  distributed  to  have  the  cross  terms 
drop  out. 

Similarly  one  obtains 


G22(w)  - (n22  + X2  + Y2)2tt6(w-wo). 


(A-8) 
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Writing  the  cross  correlation  and  ensemble  averaging  yields  the  cross 
spectrum  G^d)): 


r 2 ,(i^) 

312(w)  ” lA  ' * '+  *2e  °'  Cy  ' J 


2irfi  ((i)-uo) , 


(A-9) 


and  the  coherence  is  simply 


iu 


y - j[^j 
f12  * 


/xl  ~ M /*!  ~ x2\ 

5'  Cx  /+Y2eiwo\  Cv  / 


Jo^2  + X2  + Y2)  (n22  + X2  + Y2^1* 


This  is  Identical  to  equation  (3),  where 


(A-10) 


_ 2 2 
X2(u)  - X(f),  YZ(u)  - Y(f);  and  n^w)  - n2(u)  - N(f). 


Also,  (x^  - x2)  - d. 

When  6X  and  6y  have  deterministic  relationships,  different  results  are 
effected.  Consider,  for  example,  when  8X  - 0y  - 0,  the  observed  signals  are 

X X 

iu  t iu  (t  - -i)  iu  (t  - -^) 

Sx  - nxe  + Xe  0 c*  + Ye  y (A-ll) 


and 


X X 

id)  t id)  (t ~)  id)  (t — ) 

S2  - n2e  0 + Xe  0 cx  + Ye  ° cy 


(A-12) 


Note,  only  X and  Y are  not  independent;  however,  n^  and  n2  are  indepen- 
dent relative  to  each  other  and  to  X and  Y.  Since,  as  observed  previously, 
all  cross  terms  Involving  n^  and  n2  will  drop  out  after  ensemble  averaging, 
there  is  no  need  to  demonstrate  this  further  by  incorporating  the  intital 
phase  angles.  Also,  n^  - njCoj),  n2  * n2(d>),  X - X(id),  and  Y ■ Y(<d)  is  under- 
stood . 


The  autocorrelation  at  point 

1 is 

T/2  . 

/.  id)  T 

Tr  , *•  t 

T/i . 

(t)  - y 

/ n.n.e  0 dt  + 

a/  1 A 

f X2e  ° dt  + 

f Y2.  » 

-T/2 

-T/2 

-T/2 

lim  T -*■  « 

T/i  id)  t id)  . X1  xl\ 

id)  x id) 

/^i  _ M 

J XYe  C 

* e 0 \c  c / dt 

\ x y / 

+ f YXe  ° e ° 

\y  cx'< 

-T/2 

-T/2 

(A-13) 
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and 


G,,  (u)  * n,  2x6  (u  - u 
111  o 


Similarly, 


3 3 [-«£-?) 

) + (X  + Y )2x6(u  - uo)  + XY[e  ^x  y/ 

i-P  - iyi 

+ e \Cy  Cx/j2x6(u  - u ) 

H 


G-_(u)  - n_22x6(u  - to  ) + (X2  + Y2)2x6(u  - u ) + XY 
III  o o 


X2  V 

cx  cyJ 


+ e 


iu  ri  _ V 
vy  c*> 


2x6(u  - u ) . 
o 


The  cross  spectrum  is  readily  obtained: 


- iu 

G12(o>)  - |X  e ° 


\ cx  /+  Y2e  ^ cy  /+  XY  e °\' 


iU°*cx  " CyJ 


+ e 


iu  (-i  - — ’ 
oVc  c 

\ y x 


Also,  one  may  write 


2«5(u  - u ) . 

o 


G , (u)  •>  n.  2ti 6 (to  - u ) + 
11  1 o 
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1 xi 

-iu  ~ -iu  -1 


Xe 


-iu 


r xi  xii  r x2  vi 

» I iu  — 4iu  — -iu  — -iu  —I 

G12(u)  - [xe  ° cx  + Ye  0 Cyj  [xe  0 cx  + Ye  0 CyJ 


o cx 

+ Ye  ° 

l)  " * ~ 

-iu 

0 cx 

+ Ye  ° 

1 -iu 

^2 

(A-14) 


(A-15) 


(A-16) 


2x6  (u  - u ),  (A-17) 

o 


J_  2 

zy  2tt6  (u  - uq),  (A-18) 


2x5 (u  - u ) (A-19) 


and,  ignoring  uncorrelated  noise  sources,  as  was  done  previously  in  equation  (6), 


C12  <fa>) 


' +iu  x,/c  +iu  x, /c 

,Xe  ° 1 x + Y,  0 1 1 


)( 


-iu  x„/c  -iu  x_/c 

Xe  0 2 x + Ye  ° 2 y) 


-iu  x /c  -iux/c  |2|  -iu_x,/c 

Xe  + Ye  ° 1 y 


-iu  x,/c 

Xe  0 2 x + Ye  0 2 y 


(A-20) 
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Finally,  one  notes  that  for  a standing  wave,  cx  ■ -cv  ■ c and  equation 
(A-20)  becomes 


C12(w) 


-iu>  X_/c  4iti>_X,/c 


Xe 


+i(c  x, /c  -iw  x,/c\  / ~AUI  An  l C TAW  An 

Ue  01  +Ye  ° 1 ) Ue  0 2 + Ye  ° 2 . 

"iuoxl/c  + Ye+1“oXl/  2 Xe"iUoX2/c  + Ye+iV2/c|  2 


(A-21) 


for  X - Y, 


cos  w0(x1/c)  cos  0)o(x2/c) 

| cos  U X,  /c 
o X 


which  is  identical  to  equation  (6). 


cos  u x0/c 
o L 


(A-22) 
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